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Towards the use of asteroseismology to investigate the 
nature of dark matter 
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ABSTRACT 

The annihilation of huge quantities of captured dark matter (DM) particles inside 
low- mass stars has been shown to change some of the stellar properties, such as the 
star's effective temperature or the way the energy is transported throughout the star. 
While in the classical picture, without DM, a star of 1 M Q is expected to have a 
radiative interior during the main sequence, the same star evolving in a halo of DM 
with a density p x > 10 8 GeV cm~ 3 will develop a convective core in order to evacuate 
the energy from DM annihilation in a more efficient way. This convective core leaves 
a discontinuity in the density and sound-speed profiles that can be detected by the 
analysis of the stellar oscillations. In this paper we present an approach towards the 
use of asteroseismology to detect the signature produced by the presence of DM inside 
a star, and we propose a new methodology to infer the properties of a DM halo from 
the stellar oscillations (such as the product of the DM density and the DM particle- 
nucleon scattering cross-section). 

Key words: dark matter - asteroseismology - stars: interiors - stars: low-mass - stars: 
fundamental parameters - Galaxy: centre 



1 INTRODUCTION 

Different observations in a wide range of scales, from galac- 
tic to cosmological, suggest the existence of a new kind of 
matter, called Dark Matter (DM), formed by unknown par- 
ticles. Among the possible constituents of DM, the WIMPs, 
massive particles with non-negligible scattering cross-section 
with baryons, are considered one of the best candidates 
jBertone. Hooper fc Smj|2005h . 

Soon was realised that , if WIMPs exi s t, the y will 
accumulate inside stars (IPress fc Spergell Il985h and 
their annihilation may lead to significant changes in 
the classical picture of stellar evolution if the halo 
where the stars evolve has a very high density o f DM 
particles (iBoucmet fc Salatil 1 1983 : ISalati fc Silkl 1 19891 : 
iDearborn et al J 1 199(f )." In this context, the effects of the 
capture of WIMPs by the Sun were studied, addressing 
the prospects of helioseismology to test models that 
solved the old solar neutrino problem (|Dappen et al I I 19861 : 
iFaulkner et al.lll98(f), and to give constraints to the nature 
of DM particles dLopes, Silk fc Hansen! 120021: lLopes fc Silkl 



Recently, particular attention has been given to the first 
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ing th e prospects for their detection with the J WST tele- 
scope (|Freese et al.ll201ol ; IZackrisson et alj|201oh . Similarly, 
other authors focused on the DM effects on stars in the lo- 
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The purpose of this paper is to pave the way for the use 
of asteroseismology to provide an evidence of the footprint 
left by DM annihilation on the stellar oscillations. To do 
this, we will concentrate on solar-mass stars that evolve in 
haloes with very high DM densities, and we will show how 
asteroseismology may tell us about the properties of such 
DM haloes. 
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2 STELLAR EVOLUTION WITHIN DENSE 
DARK MATTER HALOES 

The evolution of a star within a halo of DM depends strongly 
on the ability of the gravitational field of the star to capture 
the DM particles that populate the halo. T he rate at w hich 
the DM particles are captured is given by (|Gouldlll987l l 



Cx(*) = £ / ^ 
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(1) 



where /„, (u) is the velocity distribution of the DM parti- 
cles seen by the star (which is proportional to the density 
of DM on the host halo p x and inversely proportional to 
the mass of the DM particles m x ) and Q Vi is the probabil- 
ity of a DM particle to be captured after the collision with 
an element i (which is proportional to the scattering cross- 



section of the DM particle with the nucleus i, a x 



i). The 



numerical subroutines to calculate the capture rate (equa- 
tion !]]) were adapted from the publicly available DARKSUSY 



code (jGondolo et al.ll2004T ). Our assumptions regarding this 
calculation are described in lCasanellas fc Lopes! (|2010h . 

Once DM particles are captured, they accumulate in 
a small region in the core of the star (r x ~ 0.01 R* for 
m x — 100 GeV). There, assuming that they are Majorana 
particles, they annihilate providing a new source of energy 
for the star. Capture and annihilation processes balance each 
other in a short time-scale, and consequently almost all cap- 
tured particles will be converted to energy, contributing to 
the total luminosity with L x = f x C x m Y . The factor f x , 
which in this work we assumed to b e 2/3 |Freese et al.ll20o"a ; 
llocco et al.|[200i ; lYoon et al.ll2"o"08T ) , accounts for the energy 
that escapes out of the star in the form of neutrinos. Recent 
Monte Carlo simulations suggest that the fra ction of the 
energ y lost in neutrinos may be even smaller (|Scott et alj 

l2ooa i. 

Due to this new source of energy, stars will evolve differ- 
ently from the classical picture if surrounded by a dense halo 
of DM. For very high DM densities (p x > 3 x 10 9 GeV cm" 3 
for a 1 M© star), the energy from DM annihilation pre- 
vents the gravitational collapse of the star, stopping its evo- 
lution in the pre-main-sequence phase, before the star could 
reac h enough central temperat ure to trigger hydrogen burn- 
ing (|Casanellas fc Lopesll2009r h 

For lower DM densities (10 8 GeV cm" 3 < p x < 
3 x 10 9 GeV cm -3 for a 1 M Q star), DM burning is a com- 
plementary source of energy for the star. As it is produced in 
a region much more concentrated than the nuclear burning, 
which normally extends up to 0.1 — 0.2 R*, the radiative tem- 
perature gradient (Vrad = d\nT / d\n P g ) is much steeper in 
the core of the star. Consequently, as the radiative trans- 
port is not efficient enough to evacuate all the energy in 
the central region, the star develops a convective core which 
was not present in the classical scenario without DM. The 
radius and duration of the convective core in crease when 
more energy from DM annihilation is produced (|Scott et alj 
l2009h : therefore, they depend on the density of DM in the 
place where the star evolves and on the properties of the 
DM particles. The balance between DM annihilation, nu- 
clear burning, and the gravitational energy leads to a new 
hydrostatic equilibrium with a lower central temperature. 
The star consumes its hydrogen at a lower rate, extending 
the time that it spends in the main sequence. These new 



properties allow us, as it will be shown, to provide a tool 
to infer the DM characteristics from the stellar oscillations 
using the seismological analysis. 



3 BASICS OF ASTEROSEISMOLOGY 

With the improvement on the quality of the data, asteroseis- 
mology is now becoming a precise tool to infer the properties 
of stars showing solar-like oscillations (I Michel et alj [20081 : 
iGarcfa et~a l, 2009; Bedd ing et al.ll2010T l. which are driven by 
turbulence in the superficial layers of the star. The eigenfre- 
quencies of solar-like oscillations can be approximated, for 
l/n — > (where I and n are the degree and the radial order 
of the modes), by the asymptotic expression 



(n + -+e v )u +0{v 2 ), 



(2) 



where vo = [2 ^f] _1 is the inverse of twice the time spent 
by the sound to travel between the centre and the acoustic 
surface of the star, and t v is determined by the properties 
of the surface layers. For a more in-depth explanation of the 
ba sics of the seismological a nalys is, the reader is referred 



to 



Tassoull (|l980h : iGoughl (|l985l ); lLopes fc Turck-Chiezej 



(| 19941 ). The value of vq can be estimated through the large 
separation Av n y. 



Av-n.,1 = fn,i — Vn-\,l — vq- 



(3) 



This parameter is se nsitive to t he mean density of the star: 
{M/R 3 ) 1/2 dCoxlll98Ch . while the small separation 
5v n ,l, given by 



Sfn.l — ^n,l — Vn-l,l+2, 



(4) 



is sensitive to the temperature and chemical gradient in the 
deep interior. 

In the last years, other relations between the frequencies 
of t he oscillation modes were proposed (for a re cent review, 



seelChristensen-Dalsgaard fc Houdekl (|2009h or lAerts et al.1 
<|201Ch ). broadening the diagnostic potential of seismol- 
ogy. Among the po ssible diagnostic methods of convective 
cores and envelopes ( Monteiro et al ]|l994l :l Lopes et al . 1997; 
lLopes fc GoughlboOli r we highlight the ratios between the 
small separations and t he lar ge separations developed by 
iRoxburgh fc Vorontsovl I2003T ) in order to suppress the ef- 
fects of the modelling of the near-surface layers: 

doi dio 



roi 



rio 



(5) 



where 



doi = -(fn-1,0 — 4f n -l,l + 6v n ,0 — iVn,l + fn+l,o), (6) 



dio = —-(vn-i,i - 4^n,0 



n+i,o ■ 



(7) 



The mixing of elements produced in the convective re- 
gions introduces a sharp structural variation in the border 
with the radiative regions that can be seen in the den- 
sity and sound-speed profiles. This sharp structural varia- 
tion p roduces an o scillatory signal in the frequency spec- 
trum |Gougblll99(ih whose period is related with the acous- 
tic depth of the discont inuity inside the star. Recently, 
ISilva Aguirre et alj (|2010l ) proposed the use of the ratios roi 
and rio to determine the size of a convective core by fitting 
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a sine wave to their oscillation pattern. Similarly, another 
combination of the small and large separations, 

D02 D13 



dr 



where D 



1.1+2 



(8) 

6u n ,i/(4:l + 6), was suggested by 
ICunha fc Metcalfe! |2007l i to measure the amplitude of the 
sound-speed discontinuity at the edge of a convective core. 
These seismic parameters (equations [S] and [S]) are sensitive 
to the presence of DM inside a star, given that they are 
uniquely dependent on the star's core structure and almost 
independent of the physical processes occurring in the sur- 
face layers. 



4 ASTEROSEISMIC SIGNATURE OF DARK 
MATTER PARTICLES 

To grasp the signature that the annihilation of captured DM 
particles leaves on low-mass stars we evolved a set of 1 Mq 
stars, with the same initial conditions (Z=0.018), in haloes of 
DM with different densities p x and different spin-dependent 
(SD) WIMP-nucleon cross-sections a Xt sD- Throughout our 
work, we considered fiducial values for the mass of the DM 
particles, m x = 100 GeV, and for their self-annihilation 
cross-section, < a a v >= 3 • 10 -26 cm 3 s _1 . The evolution 
of the stars was compu ted using a w ell-established stellar 
evolution code (CESAM; lMorell |l997])) used to compute so- 
phist i cated solar models for heli oseismology (|Couvidat et al.l 
120031 ; iTurck-Chieze et alj|20ic( ) and more recently used in 
the context of astero s eismic studies jKervella et al.l |2004| ; 
|Pe Ridder et alj|2006l ; ISuarez et afll201Ch . When the stars 
reached a luminosity equal to that of the Sun, a very precise 
mesh (with 1000 layers) was generated. Then, we calculated 
the frequencies of the oscillation modes of the stars using 
the ADIPLS code (|Christensen-Dalsgaardl 120081 ) . The char- 
acteristics of some of these stars are shown in Table [TJ and 
their sound-speed and density profiles, in Fig. [T] 

The accretion and the annihilation of DM particles in 
the core of the stars may change significantly their proper- 
ties. As expected, we found that the effective temperature 
of the stars that evolved in haloes with high DM densities is 
shifted to lower values (see Table [I}, due to the presence of 
a conv ective core (see Fig. [2 a), in agreement with pre vious 
works (|Fairbairn et al . 200 St lCasanellas fc Lopesll2009l ). The 



lower effective temperature and the larger radius lead to a 
decrease in the large separation Ai/„ t i (see Fig. [2]b). For a 
star with a known mass, the drop in Ai/ n: i , predicted by the 
relation Av n ,i oc A/ 1//2 i? _3//2 , is unmistakably related with 
the radius of the star. Furthermore, we also observed a drop 
on the small separation Su n<0 (see Fig. HJc), caused by a de- 
crease in the central density. The strong dependence of the 
global modes on the density profile of the star is responsible 
for that drop. 

In order to test the validity of our method, we checked 
if classical stars with similar characteristics may mimic the 
properties we described for stars evolving in DM haloes. In 
particular, we found that a star with a mass M* = 0.955 Mq 
and a metallicity Z — 0.04 reaches, near the end of the main 
sequence, the same luminosity and effective temperature as 
the star (iv) in our set (see Table [I}. At that moment, the 
radius of both stars is identical, leading to very similar great 
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Figure 1. Sound-speed (a) and density profiles (b) of 1 Mq stars 
that evolved in DM haloes with different densities p x and SD 
WIMP-nucleon cross-sections cr x sD when they reached a lumi- 
nosity L = 1 Lq (for each star, the product p x cr x is indicated in 
the legend in GeV cm -1 ). 
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Figure 2. (a) Size of the convective core, and the calculated 
scismologic parameters: (b) mean large separation (for 1=0,1,2,3), 
(c) mean small separation (for 1=0) and (d) slope of e£ro2l3, for 
1 Mq stars that evolved in DM haloes with different densities 
p x and SD WIMP-nucleon cross-sections <j Xi sd> when the stars 
reached a luminosity L = 1 Lq . 



separations (< Av nt i >= 128 pHz for star (iv) and 126 /iHz 
for the other). However, as the star that evolved without 
DM is in a later stage of evolution (X c — 0.03, while X c — 
0.38 for star (iv), the small separation, being very sensitive 
to the chemical gradient in the deep interior, allows us to 
differentiate both stars. In our case, star (iv), which evolved 
in a dense halo of DM, has a < 8v n ,o >= 7 piiz. This is 
almost double than that of the star with different M+ and 
Z (< 8v n ,o >= 4 ^iHz in that case). 
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Table 1. Characteristics of stars of 1 Mq when they reached a luminosity L = 1 Lq after 
evolving in haloes of DM with different densities p x and different SD WIMP-nuclcon cross- 
sections o x sd- The last two columns are the radius and the acoustic radius (t = J Q r — ) of 
the convective core (CC). All the stars had the same initial conditions (Z = 0.018). 





Px a x 
(GeV cm- 1 ) 


x c 


R* 
(H©) 


T eff 

(K) 


T c 
(MK) 


pc 

(g cm" 3 ) 


rcc 

(R*) 


TCC 

00 


(i) 





0.35 


1.000 


5777.5 


15.52 


148.7 


No CC 




(ii) 


10 -30 


0.38 


1.003 


5768.1 


15.66 


137.0 


0.04 


53.6 


(iii) 


10 -29 


0.40 


1.024 


5708.5 


15.25 


123.8 


0.08 


110.8 


(iv) 


2 x 10- 29 


0.38 


1.047 


5646.0 


15.78 


114.9 


0.09 


132.7 


(v) 


3 x 1Q- 29 


0.35 


1.071 


5582.4 


15.65 


108.5 


0.10 


145.6 



In addition, one of the most promising signatures of 
annihilating DM in stars is the fact that it can originate 
the formation of a convective core (unexpected in the clas- 
sical picture for stars with masses < 1.2 Mq) whose radius 
grows with the DM density p x . The convective core leaves 
a peculiar footprint in the profiles of the sound speed and 
density (see Fig. [TJ characterized by a discontinuity in the 
edge of the core. The presence of the convective core can be 
detected by the seismological analysis using a relation be- 
tween the small separation of modes with different degrees 
(and therefore with different depths of penetration inside 
the star). For that purpose low-degree modes (I = 0, 1, 2, 3) 
are chosen, because these modes are the ones that penetrate 
deep into the stellar core. 

In particular, we found that the seismological param- 
eter dro2i3 (see equation |5J is sensitive to the sound-speed 
discontinuity at the edge of the convective core and, conse- 
quently, to the characteristics of the DM halo. In Fig. 0a) 
we show the behaviour of the parameter dro2i3 for stars 
that evolve in DM haloes with different characteristics. We 
found that the absolute value of the slope of dro2i3 at high 
frequencies increases with the amplitude of the sound-speed 
discontinuity caused by the convective core, as predicted by 
Cunha and Metcalfe. Therefore, the slope of dro2i3 is di- 
rectly related with the amount of DM in the halo where the 
star evolves (see also Fig. [2]d). 

We also tested the method recently proposed by 
ISilva Aguirre et all |2010l) designed to estimate the size of a 
convective core in 1.5 Mq stars. We found that the period of 
the sinusoidal fit to the ratios roi and rio (see Fig. 0b) does 
not match exactly the acoustic radius of the convective cores 
(see Table [TJ, most probably because we are applying this 
method to stars of mass 1 Mq. However, the ratios roi and 
no have a great sensitivity to the amplitude of the sharp 
variation of the sound speed caused by the annihilation of 
DM particles inside the star. We conclude that these ratios 
may be used in the future as a stellar probe to confirm the 
presence of DM in the neighbourhood of low-mass stars. 

If enough radial modes are identified with the preci- 
sion presently achieved by space-based telescopes as CoRoT 
(a relative error on t he individual frequencies of ~ lO -4 
|Deheuvels et al.ll2010t )), then our method will allow the dis- 
crimination between haloes of DM with different character- 
istics. To illustrate this point, we plotted in Fig. 3 the error 
bars on dro2i3, J"oi an d r io for star (iii) derived from the 
mentioned uncertainty (I0~ 4 v) on the determina tion of the 
frequencies, as done bv lCunha fc Metcalfe! (|2007l ). 
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Figure 3. (a) The seismological parameter dro2i3 and (b) the 
ratios roi and rio for stars that evolved in DM haloes with dif- 
ferent densities p x and different SD WIMP-nucleon cross-sections 
a x-SD when they reached a luminosity L = 1 Lq (for each star, 
the product p x cr x is indicated in the legend in GeV cm -1 ). Er- 
ror bars are shown for star (iii) assuming a relative error on the 
identification of the frequencies of 10" 4 . 



5 DISCUSSION AND CONCLUSIONS 

In this paper, we have presented a new methodology to- 
wards the use of asteroseismology to prove the presence 
of DM in the location where a star evolves. For a main- 
sequence star of 1 Mq evolving in a DM halo with a density 
p x > 10 s GeV cm" 3 (assuming a XtS D = lO" 38 cm 2 ), the 
annihilation of captured DM particles on its interior leads 
to decreases in the large and small separations, when com- 
pared with the same star in the classical scenario without 
DM, which are related to changes in the global properties 
of the star. Furthermore, the highly concentrated produc- 
tion of energy by DM annihilation creates a convective core 
which is not present in the classical picture. This convective 
core leaves a discontinuity signature in the sound-speed and 
density profiles which can be detected by the analysis of the 
stellar oscillations. 

We have shown that seismological parameters such as 
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Figure 4. DM densities at which 1 Mq stars with a luminos- 
ity 1 L0 are expected to show strong signatures on the seis- 
mological parameters Av, 8u, dro2i3, rgi an d r l0 (see text) due 
to the annihilation of DM particles with different characteristics 
(m x ,o x s£)) in their interior . In the particular case of the Galactic 
Centre (GC), the DM densities in the Figure (from top to bottom) 
are expected at a distance from the GC of 0.1 pc, 0.04 pc and 
0.02 pc. following the adiabatically contracted profile of Bcrtone 
& Merritt. The grey lines are the present limits from direct detec- 
tion experiments: XENON10 (dotted), PICASSO (dashed) and 
COUPP (solid), and the grey region is the DAMA/LIBRA al- 
lowed region. 



d)"02i3 and the ratios roi and rio are very sensitive to the size 
of the convective core, which is determined by the density of 
DM p x where the star evolved and by the scattering cross- 
section of the DM particles off nuclei a x . Consequently, this 
relationship may be used in the future to help in the deter- 
mination of these parameters (or at least to their product, 
p x o~ x ) and to provide a stellar probe that identifies the pres- 
ence of self-annihilating DM. 

The method presented in this paper is valid for haloes 
with very high DM densities. In Fig. 0] we show the DM 
densities at which a 1 Mq star with a luminosity 1 Lq 
is expected to have a small separation 25% smaller than 
that in the classical picture, because of the annihilation of 
DM particles with different characteristics (m x ,a x - +Pt sD) 
in its interior. If found, this kind of star will have also 
strong signatures on the seismic parameters Av, dro2i3, ro± 
and rio when compared with a star with the same lumi- 
nosity that evolved without DM. In the same Figure are 
also shown the current limits on a^ ^ v ,sD from the direct 
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detection experiments XENON10 I Angle et all l200&t), PI 



CAS SO (|Arehambault et al.ll2009l ). COUPP |Behnke etal] 
120081 ') and the allowed region from the DAMA/LIBRA ex- 
periment |Savage et al.ll2009l ). 

The extreme DM densities shown in Figure [4] may be 
present within the inner parsec of our Galaxy, according 
to models that account for the effect of the baryons on 
the DM halo via adiaba tic contraction jBlumenthal et all 
Il986l ; iGnedin et al.|[20ol ). For i nstance, following t he ad i- 
abatically contracted profile of iBertone fc Merrittl (|2005h , 
DM densities as high as p x = 10 8 GeV cm -3 are ex- 



pected at 0.1 pc from the Galactic Centre (GC). Even higher 
DM densities may be present at the GC if a hypothetical 
spike is formed due to th e influence of the central black 
hole (|Gondolo fc Silklll999l ). However, as other models pre- 
dict lower cen tral DM densities (the so-called core models 
|BurkertJI 19951 )1. the open questions about the DM halo pro- 
file at the inner parsec of our Galaxy are still far from being 
solv ed (for a r ecent review on this topic, see iMerrittj (|201oh 
and Ide Blokl (|2010T l). In this sense, the method proposed 
here may provide a complementary tool to help in the dis- 
crimination of different models. Other possible locations of 
environments with such high DM densities are the dwarf 
spheroidal galaxies around the Milky Way (|Dekel fc Silkl 
1 19861 ; iKormendv fc Freemanll2004l ; iDiemand et al.ll2007l ). ~ 



The precision required for our analysis is similar to the 
one achieved by present asteroseismic missions in observa- 
tions of one hundred days. Nevertheless, the most likely 
place to find the kind of stars described here is near the 
centre of our Galaxy, where the distance and the presence 
of dust makes the observations difficult. These difficulties 
encourage us to extend our study to more massive and lu- 
minous stars in a future work. Future technical improve- 
ments in the observations of the GC and of the Milky Way 
dwarf galaxies may open the possibility of using the method 
proposed here to investigate the nature of DM. 
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